Coherence properties of nanofiber-trapped cesium atoms 
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We experimentally study the ground state coherence properties of cesium atoms in a nanofiber- 
based two-color dipole trap, localized ~ 200 nm away from the fiber surface. Using microwave 
radiation to coherently drive the clock transition, we record Ramsey fringes as well as spin echo 
signals and infer a reversible dephasing time of T2 = 0.6 ms and an irreversible dephasing time of 
T2 = 3.7 ms. By modeling the signals, we find that, for our experimental parameters, T2 and T2 are 
limited by the finite initial temperature of the atomic ensemble and the heating rate, respectively. 
Our results represent a fundamental step towards establishing nanofiber-based traps for cold atoms 
as a building block in an optical fiber quantum network. 
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Over the past years, hybrid quantum systems have at- 
tracted considerable attention [T]. In the specific case 
of light -matter quantum interfaces [IHlj, they combine 
the advantages of photons for transmitting quantum in- 
formation and of long-coherence-time systems, such as 
dopant ions in crystals, NV centers, quantum dots, sin- 
gle trapped neutral atoms and ions, and atomic en- 
sembles, for storing and processing quantum informa- 
tion and for realizing long-distance quantum communi- 
cation [5j. In the context of quantum networks [6], it 
would be highly desirable to connect these matter-based 
storage and processing units via optical fiber links. A 
promising approach towards the realization of such fiber- 
based quantum interfaces consists in coupling cold neu- 
tral atoms to photonic crystal fibers [3{9] . Another tech- 
nique with high potential involves trapping and inter- 
facing cold atoms in the evanescent field surrounding 
optical nanofibers. Using the optical dipole force ex- 
erted by a blue- and a red-detuned nanofiber-guided light 
field [lOl [11] , two-color traps have been demonstrated ex- 
perimentally with laser-cooled cesium atoms [TU |T3] . 

In order to implement quantum protocols with atoms 
coupled to nanophotonic devices, good coherence prop- 
erties are a prerequisite but cannot be taken for granted: 
various effects, like Johnson noise [13] or patch po- 
tentials [15] may occur and hamper long coherence 
times [H]. When coupling to optical near-fields, this is 
all the more critical because of the small atom-surface 
distance of typically a few hundred nanometers. Here, 
using Ramsey interferometry as well as spin-echo tech- 
niques, we measure, to the best of our knowledge for 
the first time, the reversible and irreversible dephasing 
times of atoms in such an environment. Specifically, we 
experimentally characterize and model the ground state 
coherence of the clock transition of cesium atoms stored 
in the nanofiber-based two-color trap realized in 12J. Re- 
markably, the inferred coherence times extend up to mil- 
liseconds. 

The experimental set-up is sketched in Fig. [T^ and is 
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FIG. 1. a. Sketch of the experimental set-up including the 
tapered optical fiber, the trapping, probe and push-out laser 
fields, the microwave antenna, and the single-photon counter 
(SPCM). b. Orientation of the plane of the quasi-linear po- 
larizations of the blue- and red-detuned trapping fields, the 
atoms, and the magnetic offset field in a plane transversal to 
the fiber axis. 



described in detail in [T?, Tf^ . Cesium atoms are trapped 
in the evanescent field surrounding the nanofiber-waist 
of a tapered optical fiber. The atoms are located 
about 200 nm above the nanofiber surface in two dia- 
metric one-dimensional arrays of potential wells, with 
at most one atom per trapping site. Localization in 
the three (radial, azimuthal and axial) directions is 
achieved with trap frequencies of (200, 140, 315) kHz. In 
order to drive transitions between the hyperfine ground 
states of the trapped atoms, we use a tunable microwave 
(MW) field at a frequency of 9.2 GHz. In the follow- 
ing, we limit our study to the so-called clock transi- 
tion between the states |e) = |65i/2, F = 4, mp = O) and 
\g) = \QSi/2iF = 3, rap = O). This \g) — >■ |e) transition 
only exhibits a quadratic Zeeman shift and can be se- 
lectively addressed by applying a homogeneous magnetic 
offset field. Bog. The latter is oriented perpendicular to 
the plane containing the atoms, see Fig. [TJa. 

We load the cesium atoms from a magneto-optical trap 
into the nanofiber trap via an optical molasses stage [H] . 
In order to then prepare an ensemble of atoms in \g), we 
perform a state purification sequence: first, all atoms 
are optically pumped into the F = 4, hyperfine ground 
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FIG. 2. a. Rabi oscillations: Probability pe(fp) for an atom 
to be in state |e) after a resonant MW pulse of length tp. The 
data points correspond to an average over 80 experimental 
realizations (scaling factor a = 0.015). The solid line is the 
result of a global fitting routine (see text), b. Ramsey in- 
terferometry: Data averaging and fitting procedure like in a.; 
the probability p'^'^'^\t) is plotted as function of the time 
interval r between the n/2 pulses (scaling factor a — 0.018). 
c. Vibrational state-dependent differential light-shift (5is(n) 
(squares) and histogram of the phonon number probability 
distribution P{n, T) calculated at T = 71 (xK. 



state, then the atoms in |e) are selectively transferred to 
\g) by a MW pulse, and, finally, the remaining atoms in 
F = A are removed from the trap by means of a cr+- 
polarized "push-out" laser beam resonant with the light- 
shifted = 4 — )■ F' = 5 transition of the cesium D2 
line p[H]. We characterized the efficiency of this purifica- 
tion procedure and found that at least 62 % of the re- 
maining atoms end up in state \g) while the population in 
the F — A manifold is negligible. In our nanofiber-trap, 
spin flips from \g) to adjacent Zeeman substates may oc- 
cur because the atoms move in a light-induced fictitious 
magnetic field that exhibits a strong gradient of several 
Gauss per micrometer [191 [20]. We observed, however, 
that these spin flips are suppressed for a magnetic offset 
field of Boff > 3 G. Therefore, a field of 3 G is applied for 
all coherent manipulations of the atomic spins. 

We first record Rabi oscillations between states \g) and 
|e) by applying a resonant MW pulse of variable dura- 
tion. The transfer probability — N^/Ntot is given by 
the ratio of the number of atoms that were transferred 
to |e), denoted N^, and of the sum of the numbers of 



atoms in |e) and \g), denoted A^tot- Given the fact that 
the F = A manifold is depleted after the purification and 
that the MW selectively drives the \g) — > |e) transition, 
we identify with the number of atoms in = 4 and 
measure the latter by recording the transmission spec- 
trum of a quasi-linearly polarized fiber-guided light field 
scanned over the light-shifted F = A ^ F' = 5 transition 
of the D2 line [12] after setting Bos to zero. Because of 
shot-to-shot fluctuations of N^ot: the latter has to be de- 
termined for each experimental run. However, we do not 
have direct access to this quantity because the probing 
on the F = A F' — 5 transition leads both to loss 
of atoms from the trap and to transfer of atoms from 
|e) to \g) by Raman scattering. Instead, we measure the 
number of atoms in the trap prior to the state purifi- 
cation procedure, Moad, which is proportional to iVtot- 
The ratio a = iVtot/Moad is found to be stable and is 
determined from the analysis presented below. 

The inferred probability = iVe/(Q;A'ioad) for an atom 
initially in state \g) to be in |e) after the MW pulse is 
shown in Fig. [2^. The observed Rabi oscillations (reso- 
nant Rabi frequency r2o/27r ~ 17 kHz) exhibit no visible 
damping over 200 \xs. Longer measurements (not shown) 
reveal an exponential reduction of the amplitude of the 
oscillations which decay to Pe — 0.5 in a characteristic 
time of (3.4 ± 0.2) ms. This damping originates from 
fluctuations of the power of the MW source. 

In order to study the coherence properties of the 
atomic ensemble during free evolution, Ramsey interfer- 
ometry is performed: Two successive 7r/2-pulses, sepa- 
rated by a time r, are applied and the transfer probabil- 
ity pi^™''(r) is measured, see Fig. The frequency of 
the MW source is detuned by —5 kHz with respect to the 
frequency used to record the Rabi oscillations. The ob- 
served Ramsey fringes reveal a coherence time of several 
100 |J.s. Their decay can be explained by inhomogeneous 
broadening of the clock transition which is caused by 
the flnite temperature of the atomic ensemble in combi- 
nation with the dependence of the transition frequency 
on the vibrational state of the atoms due to differential 
light-shifts [H]. These differential light-shifts are pro- 
portional to the intensity of the trapping fields and thus 
depend on the position r of the atom in the trap. We 
denote the frequency of the clock transition at position 
r as uj{r) and the corresponding differential light-shift as 
Sis{r) = aj(r) — Wo, where ujq is the transition frequency in 
free space. The variation of Sis{r) in the azimuthal and 
axial trapping directions is significantly smaller than in 
the radial direction. Therefore, in the following, only the 
radial dependence is taken into account. For the calcula- 
tion of the motional eigenf unctions, we make the simpli- 
fying assumption that the motions in the three directions 
are decoupled and take the strong anharmonicity in the 
radial direction into account [TT]. An atom in a certain 
radial vibrational state \n) will experience the differen- 
tial light-shift d\s{n) = (n| 5\s{f) \n), see Fig. [2];, where n 
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is the phonon number and the trap parameters are the 
same as in [T2] . 

Thus, the density matrix of the atomic pseudo spin-1 /2 
just before probing, calculated in the MW field rotating 
frame, depends on n and reads 

=XR(t;n)poXR^(i;ri), (1) 
A^R(t; n) = 7^^/2 UQ{t; n) 7^^/2, 
Un(t\n) = exp [-i{6\s{n) - (5mw)^ 0-^/2] , 

where is the Pauli matrix, po — Iff) (.91? ^-nd (5mw = 
■^MW ~ '^o is the detuning of the MW field with respect 
to the clock transition in free space. Here, A^R(t; n) is the 
evolution operator corresponding to the total experimen- 
tal sequence, while Z^o(t; n) is the free evolution operator. 
The operator TZq corresponds to a rotation of the atomic 
spin by the angle about the y-axis of the Bloch sphere. 
We assume a Boltzmann distribution for the initial occu- 
pation probability T) of the vibrational states which 
is plotted in Fig. [2]; for a temperature of T = 71 |a.K. For 
a given P{n,T), the probability pi^™''(T) of detecting 
an atom in |e) for a time delay t then reads 
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In Fig. [Id, we plot pi^"""^ (r) for T = 71 and find 
that our model describes the experimental data well. We 
define the reversible dephasing time as the time at 
which the amplitude of the Ramsey fringes is reduced by 
50%, yielding = 0.6 ms. 

In order to reverse the dephasing caused by the inho- 
mogeneous broadening of the transition, the well-known 
spin echo (SE) technique is employed: A 7r-pulse is ap- 
plied in between the two 7r/2-pulses at a time tocho/2 
after the first 7r/2-pulse. In Fig. |3j we plot the mea- 
sured probability pi^^'*(T; iocho) of detecting an atom in 
|e). Spin echoes are apparent even for iocho ^ 

The amplitude of the echo signal in Fig. [3] decreases 
within a characteristic time of a few milliseconds, indicat- 
ing the presence of an irreversible dephasing mechanism. 
A general way to model the latter relies on the multi- 
plication of the coherences (off-diagonal elements) of the 
density matrix p by a time-dependent factor < C(t) < 1 
(C(0) = 1). This is done by introducing the superopera- 
tor C{t) defined by 

for(..)e{e,5r. (3) 

The function C{t) thus provides a measure of the co- 
herence of the superposition [211 [22|. The atomic spin 
density matrix at the end of the SE sequence reads 
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FIG. 3. Probability pI^^\t) for an atom to be in |e) after a 
SE sequence. The data points correspond to an average over 
80 experimental realizations (scaling factor a = 0.018). The 
solid lines are the result of the global fitting routine. The 
legend in each panel indicates tccho. 



The evolution operator A^r in Eq. ([ij has been replaced 
by the evolution superoperator AdsE- Moreover, C{t) 
is assumed to commute with the operators Uq and TZt^, 
which is fulfilled for ideal 7r-pulses and a good approxi- 
mation in the case of small imperfections. For simplicity, 
we take C (t) to be constant within the measurement time 
intervals |r — ^ochol ^ 0.5ms and, thus, replace C(t) in 
Eq.(|4| by C(techo): 



(5) 
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From this 
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expression, the transfer 

pr"'(T; tccho) is calculated as in Eq. (H. 

The solid lines in Fig. [2^, b, and in Fig. [3] are the re- 
sult of a global fit of our model to the Ramsey fringes as 
well as to the SE signals. Rabi oscillations recorded over 
~ 5 ms are also included in the global fit. The fitting pro- 
cedure disregards vibrational levels beyond nmax — 70, 
which is well-justified for typical temperatures of the en- 
semble in the nanofiber-based trap. Imperfections of the 
TT- and 7r/2-pulses due to the dependence of the rota- 
tion angle on n via Sis{n) are included. Furthermore, a 
phenomenological phase of the second 7r/2-pulse that ac- 
counts for a drift of the mean clock transition frequency 
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during the sequence due to heating (see below) is added 
as a fit parameter. The other fit parameters are the initial 
temperature Tq, the MW detuning ^mw in the Ramsey 
and SE experiments, the resonant Rabi frequency Hq, the 
coefficients C(techo), and the ratio a — iVtot/Moad- The 
experimental data and the theory shown in Figs. [2] and [3] 
agree very well. Moreover, the obtained initial tempera- 
ture of To = (71 ± 4) \xK is in reasonable agreement with 
what is expected from independent measurements of tem- 
perature and heating rate |T7J[53]. We note that the error 
bar of Tq is only statistical and does not include possible 
systematic deviations due to, e.g., the finite accuracy of 
the calculation of the differential light shifts. The qual- 
ity of the global fit together with the realistic range of 
the obtained fit parameters lead us to conclude that the 
measured reversible dephasing time is, indeed, a re- 
sult of the joint effect of the finite initial temperature of 
the trapped atoms and the position-dependent differen- 
tial light-shift of the clock transition. 

We now model the decay of the spin echo amplitude, 
quantified by the fitted parameters C(tocho), plotted in 
Fig. |4] It has been shown previously that the heating rate 
of atoms in our nanofiber-based trap currently lies in the 
mK/s-range [23]. We infer the effect of heating on the 
coherence of the trapped atoms and the resulting transfer 
probabilities p'^^\t; techo) by calculating the time evolu- 
tion of the density matrix p of the joint system, consist- 
ing of the atomic spin and the vibrational phonons. We 
numerically solve the master equation dp/di = Cp writ- 
ten in the Lindblad form [21], with the initial state p{0) 
corresponding to the tensor product of po = \g){g\ and 
the thermal state X^n^l^-'^o) I"-) ("-I- The Liouvillian C 
is given by 



+ K L^pLl - ^ {lIl^p + pLlL^ 



(6) 



where i/o = (<^is("-) ~ ^mw)<7z/2 is the free evolution 
Hamiltonian, the operator L_ = a (i+ = a^) is the 
phonon annihilation (creation) operator, and fi = a) a 
is the phonon number operator. This Liouvillian de- 
scribes the process of the atoms being heated due to a 
jitter of the center of the trap [25^. From the resulting 
Pe^^' (t; iccho), we deduce the expected time evolution of 
the coherences Ch(iccho)- Here, in view of the typical 
temperatures of the atomic ensemble and in order to 
reduce computation time, the phonon Hilbert space is 
truncated at n,nax — 30. Figure |4] shows the results ob- 
tained for K — 350 s^"'^, corresponding to a heating rate of 
7 ~ 3 mK/s which is close to the one reported in [23], and 
To = 71|J.K, taken from the global fit of the data in Figs. [2] 
and [3] The theoretical prediction matches very well with 
values of C(iocho) extracted from the experimental data. 
We therefore conclude that for our current experimen- 
tal parameters, the irreversible dephasing is dominated 
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FIG. 4. Irreversible decay of the atomic ground state coher- 
ence. Black squares: fitted parameters C(techo)- The error 
bars are obtained by the standard error propagation method. 
Red dots: numerical simulation results obtained with a heat- 
ing rate 7 ~ 3 mK/s and an initial temperature To = 71 (iK. 



by heating of the atoms. We define the irreversible de- 
phasing time T2 as the half width at half maximum of 
C'h(tocho) and obtain T2 = 3.7ms. 

Summarizing, we showed that ground state coherence 
times in the range of milliseconds can be achieved in our 
nanofiber-based atom trap |12j with atoms at a distance 
of only ~ 200 nm from the hot silica fiber surface [26j . 
Currently, the coherence times are limited by the finite 
initial temperature of the atomic ensemble and the heat- 
ing rate. Both can in principle be further decreased by 
technical measures like Raman side-band cooling [27 l l28 j . 
This constitutes a decisive result towards establishing 
nanofiber-based quantum interfaces as practical building 
blocks in an optical fiber quantum network. In partic- 
ular, they pave the way towards the realization of fully 
fiber-integrated quantum memories which could also be 
operated as highly efficient photon number-resolving de- 
tectors [23 . 
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